A mathematical model is developed for the SOLUTION POLYMERIZATION of methyl methacrylate (MMA) in a batch reactor. The model includes the whole process, In POLYMERIZATION part, model includes chain transfer and termination by both combination and disproportionation, Variation in reactant mixture properties is included in the model, The glass and gel effects are included by CCS (Chiu, Carratt and Soong) model, Energy balances are applied to the reactor, oil bath, and temperature sensors. Therefore, the model can attribute to the non-isothermal process. The final model consists of fourteen nonlinear stiff differential equations. Output of the model is monomer conversion, initiation concentration, radical concentration, dead polymer concentration, volume, and reactor and jacket temperatures. Agreement of model prediction and experimental results is very good. This model can be used for optimization and control studies because it includes the whole process.